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Abstract 

We analyze a macroscopic model with a maximal density constraint 
which describes short range repulsion in biological systems. This system 
aims at modeling finite-size particles which cannot overlap and repel each 
other when they are too close. The parts of the fluid where the maximal den- 
sity is reached behave like incompressible fluids while lower density regions 
are compressible. This paper investigates the transition between the com- 
pressible and incompressible regions. To capture this transition, we study 
a one-dimensional Riemann problem and introduce a perturbation problem 
which regularizes the compressible-incompressible transition. Specific dif- 
ficulties related to the non-conservativity of the problem are discussed. 

Keywords: Congestion, Riemann problem, incompressible-compressible tran- 
sition, clusters dynamics, gregariousness, steric constraints 
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1 Introduction 



We consider a macroscopic model of self-driven particles which describes the 
dynamics of a large number of social interactive agents. More specifically, we 
are interested in modeling short range repulsion effects due to the fact that finite- 
size agents (e.g. sheep in a herd) cannot overlap (non-overlapping or steric con- 
straints). To this aim, we derive a hyperbolic problem with a density constraint as 
a limit of an unconstrained system with a repulsive force which turns on suddenly 
when the density becomes close to the maximal one. The limit model requires 
transmission conditions at the transition between an unclustered region (where 
the maximal density is not reached) and a clustered region. In unclustered regions, 
the fluid is compressible while it becomes incompressible in the clustered ones. 
Therefore, this paper aims at providing a description of this transition between a 
compressible and an incompressible fluid. Unfortunately, the formal perturbative 
approach which we implement does not directly provide information about these 
transmission conditions. In order to retrieve this information, we rigorously ana- 
lyze special solutions of the perturbation problem: the Riemann problem. These 
solutions are explicitely known and allow us to carry out the limit rigorously and to 
recover the required transmission conditions. We postulate that these conditions, 
which are rigorously proven only for Riemann problem solutions, do extend to all 
solutions. However, being non-rigorous for general solutions, these conditions are 
stated as formal conditions in the "formal statement 1" below, which constitutes 
the main result of the present paper. Still, the rigorous analysis of Riemann prob- 
lem solutions is quite technical and the proofs of many statements are deferred to 
appendices. 

The modeling of biological systems undergoing flocking or herding dynam- 
ics has been the subject of a vast literature. A first class of models relies on the 
alignement interaction between neighbouring self-propelled particles. The sim- 
plest of these models is an individual-based (or microscopic) model proposed by 
Vicsek [25,44]. A macroscopic version of the Vicsek model is derived in [21] and 
a coUisional Vicsek model is proposed in [8]. A variant of the Vicsek model has 
been proposed by Cucker and Smale [17, 18] (see also [12,26] for recent mathe- 
matical results). By incorporating long-range attractive and short-range repulsive 
forces to the Vicsek model, one obtains the three zones model of Aoki [1,16, 39], 
originally devised to describe fish schools. Models with repulsive-attractive in- 
teraction only (without alignement interaction) have been studied in [14, 23, 32]. 
Such models have been used for pedestrian interactions [28,34]. Other kinds of 
macroscopic models of drift-diffusion type have been analyzed in [27, 31,41,43] 
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and different hyperbolic models are compared in [5]. For biological reviews, we 
can refer to [15,33]. 

As outlined above, we focus here on the congestion constraint: animals or 
individuals cannot overlap (steric constraint). As a consequence this congestion 
constraint leads to the existence of a maximal density p*, which cannot be ex- 
ceeded inside the flock. This problem has been analyzed before and schemati- 
cally two methods have been proposed. A first one consists in modeling repulsion 
through forces or diffusion terms [14,23,31,32,43]. However, in this approach, the 
individuals are point particles and their finite size is not explicitely described. So 
the maximal density constraint is not explicitely taken into account. To explicitely 
take this maximal density constraint inco account, in [29,30], the authors have de- 
velopped an alternative approach: the particles are first evolved freely over one 
time step and then projected towards the "closest" admissible non-overlapping 
configuration. This leads to non-local interaction between the particles which 
contradicts the local character of the interactions in most biological systems. By 
contrast, we developped a third route inspired by multi-phase flows [10] and traf- 
fic jam modeling [6,7]. The repulsive force is modeled by a nonlinear pressure 
law p{p) which becomes singular as the density approaches the maximal density 
p*. Additionnally a small parameter e allows to describe the fact that the regu- 
larized pressure is very small of order £ as long as the density p is smaller than 
p* and turns on suddenly to a finite or even large value when p becomes close 
to p*. In the limit e ^ of this model, two distinct phases appear: a pressure- 
less compressible phase which describes free motion in unclustered regions and 
an incompressible phase which describes the motion inside the clusters. The ma- 
jor difficulty is to find the transmission conditions between the compressible and 
incompressible phases. 

The present paper is a multi-dimensional extension of the methodology pre- 
sented in [6,7, 10] for multi-phase flows or traffic. However, an additional diffi- 
culty arises due to the non-conservative character of the original hyperbolic model. 
Indeed, momentum is not a conserved quantity because the particles in the under- 
lying particle system are self-propelled particles which have constant (in-time) 
and uniform (in-space) velocities. Therefore, the model which is at the starting 
point of this paper is a non-conservative hyperbolic system which as such presents 
an ambiguity in the definition of weak solutions. We will show that this ambiguity 
can be partly removed for one-dimensional Riemann problem solutions. We be- 
lieve that the strategy developped in this paper to analyze congestion effects can 
apply to other systems such as bacteria populations [35], economic systems like 
supply chains [2] or physical systems like granular materials [4, 38]. 
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The organization of this article is as follows. In section 2, we present the 
perturbation model and its limit. We also provide the connection conditions be- 
tween the compressible and incompressible phases of the limit model, which are 
the main result of the paper. A remark on collision of clusters is also formulated. 
With these informations, we show that the available information is sufficient to 
provide a well-defined dynamics at least in the case of a single cluster. Section 3 
is devoted to the study of the one-dimensional Riemann problem for the perturbed 
problem and the limits of its solutions as e — 0. As stated above, this analysis 
provides a strong support for (but not a proof of) the postulated transmission con- 
ditions at the compressible-incompressible interface which are provided in section 
2. Appendix A provides a formal derivation of the initial model from an individ- 
ual based model with long-range attractive and short-range repulsive interactions, 
which describes the aggregation of gregarious animals like sheep. Appendices B 
to E provide proofs of technical lemmas and theorems needed in the analysis of 
the Riemann problem. 

2 Model and goals 

2.1 The model and its rescaled form 

Our starting point is the following model, written in dimensionless form: 



where p = p{x,t) is the particle density and Q. = Q{xj) is the particle velocity. 
The problem is posed on the 2-dimensional plane x G and ? > is the time. 
The velocity ^(xj) E M? is supposed to satisfy the normalization constraint 



Therefore, Q.{x,t) E S , the unit sphere, at any point in space-time. The function 
p{p) is an increasing function such that p{p) ~ when 7^ 1 and p{p) +00 
when p ^ p* where p* is the so-called congestion density. In this paper, we will 
consider 



d,p + VriP^)=0, 

dta+{a-v^)a + {id-si®Q.)Vp{p) = o. 



(2.1) 
(2.2) 



nix,t)\ = 1, V3ceM^ \/t>o. 



(2.3) 
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(2.4) 



4 



for simplicity but any other function with similar behaviour would lead to sim- 
ilar results. The operators V^- and {Q. ■ V^) are defined, for a vector field A = 
(Ai,A2)(x),by 

V^.A = d,^Ai+d,^A2, (2.5) 

{a-v^)A = {{aA, +a2d,^)Au{nid,^ +a2d,-M2V, (2.6) 

where T denotes the transpose operator. Finally, (Id — Q ® ^) is the projection 
matrix onto the line spanned by Q.^, where Q.^ is the vector Q. rotated by the angle 
7r/2. Alternatively, we have, for a vector A: 

{ld-Q^a)A=A-{a-A)A, (2.7) 

where {D. -A) is the dot product Q. A = Q.iAi +Q.2A2. 

We show in appendix A that this model well describes the behaviour of a sys- 
tem of particles subjected to long-range attraction and short-range repulsion in 
the spirit of a model proposed by Aoki [1] or Couzin et al [16] for modelling gre- 
gariousness and swarming. More precisely, in appendix A, we derive this system 
from such a particle system through successive changes of scales via mean-field 
and hydrodynamic theories. In the form (12.21) . we have dropped the force term 
describing long-range attraction. Indeed, this force term would add the quantity 

(Id - 12 (g) Q.) la at the right-hand side of (IZ2l) . with 

E* lKa{\y-x\){y-x)p{yj)dy 
^"^'^ lKa{\y-x\)p{yj)dy ' 

where Ka is a bounded positive kernel. This terms does not add any differential 
operator and all the subsequent analysis will stay unaltered by adding this term. 

Our main concern is the study of the congestion effects brought by the singu- 
larity of p{p) near the congestion density p*. Indeed, a herd of animals can be 
viewed, at large scales, as a domain of space where the density p is close to the 
saturation density p*. Therefore, the geometrical domain occupied by the herd at 
time t can be identified to a set H; = {x eM?\p* - 5p < p{x,t) < p*} where the 
parameter 5p > must be suitably tuned. Therefore, with the initial model (12.11) . 
(12.21) . the definition of a herd depends on an arbitrary parameter Sp, which makes 
it ambiguous. 

A way to unambiguously define the herd is to force the system (I2.1I) - (I2.2I) to 
make clear-cut phase transitions from unclustered p < p* to clustered p = p* 
phases. In the spirit of the works [6,7,20] for traffic, this can be achieved in an 
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(a) e = 1 (b) £ = IQ-^ 



Figure 1: The "potential" for repulsive interaction p{p) (left) and £p{p) (right) 
after scaling by a small parameter £ = 10^^, with 7=2 and p* = 1. From the 
right picture, it is clear that the repulsive interaction turns on only when p is very 
close to p* 



asymptotic regime which amounts to supposing that there is merely no repulsive 
interactions at all as long as p < p*, and that repulsive "pressure" forces turn on 
suddenly when p hits the congestion density p*. This can be done by rescaling 
p{p) into £p{p) where £ ^ 1 is a small parameter. In this way, repulsive in- 
teractions are 0(e) as long as p < p*, but become 0(1) when p = p* (see fig. 

m. 

Biologically, this assumption amounts to saying that the animals do not change 
their directed motion by the presence of their neighbours unless they touch them 
and need to modify their trajectory to bypass them. The parameter £ <C 1 is re- 
lated to the time scale at which this change of trajectory occurs and is therefore 
supposed small. Let us also note that our model considers that all animals move 
with speed unity and never stop. Obviously the model will require improvements 
by taking into account the fact that a certain fraction of animals are steady, while 
foraging or resting. 

Therefore, our main concern in this paper is the study of the following pertur- 
bation problem: 

drp' + Vr{p'^')=0, (2.8) 
dtO.' + {a'' ■Vj.)Q.' + e{ld-Q.' (^Q.')V^p{p') =0, (2.9) 

= 1. (2.10) 

We will be interested in the formal limit £ ^ 0. A rigorous theory of this type 
of problems is unfortunately still out of reach up to our knowledge. In the fol- 
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lowing section, we show that the limit £ ^ leads to a phase transition between 
compressible and incompressible regimes. 

2.2 The singular limit e ^ 0: transition between compressible 
and incompressible motion 

As £ — i> 0, £p{p^) becomes significant only where the convergence p*^ ^ p* is fast 
enough. Therefore, in the limit, either ^ p < p* and £p{p^) ^ or p^ ^ p* 
and £p{p^) p with p possibly non zero. In other words, the equation (p* — 
p)p = holds in the limit. If additionally p < +oo, straighforward inspection 
shows that 



Therefore, the formal limit £ ^ of system (|2.8|) - (|2.9I) - (I2.10I) is given by the 
following system: 



In the non-congested domain p < p*, the system reduces to a pressureless com- 
pressible gaz dynamics model with a speed constraint 



This system describes the behaviour of the system outside the congested region. 
It is a compressible system. Biologically, it describes the behaviour of dispersed 
animals outside the herd. Mathematical studies of this system are outside the 
scope of this article and the reader can refer to [9] for standard pressureless gas 
dynamics models (without speed constraint). We note that this system exhibits 
vacuum regions where p = as it will be seen below. 

2.3 Study of the congested region 

The congested part of the flow is defined as the region where the congestion con- 
straint p = p* is reached. Biologically, it defines the domain of space occupied by 



p*-p^ = 0{£r). 



(2.11) 



dtp+yr{p^)=0. 

dtSi+a-v^a+{id-a0a)v^p = o, 
\a\ = 1, 

(p*-p)p- = 0. 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



dtP + VrP^ = 0, 

\a\ = 1. 



(2.16) 
(2.17) 
(2.18) 
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the herd. Its connected components will be called "clusters". In the congested do- 
main, system ( |2.12|) -( |2.15l) turns into an incompressible Euler model with speed 
constraint: 



dtQ+n- v^a + {id-Q® n)ysp = 0. 
1^21 = 1, 

p = p\ 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



We first note that smooth incompressible vector fields of constant norm in M? 
have a very special structure which is outlined in the following. 

Proposition 2.1 Let Q.(x) be a smooth vector field on a domain C with 
values in and which satisfies the incompressibility constraint V^ Q. = Q. Then 
the integral lines of Q.^ are straight lines and Q. is constant along these lines 
( where is rotated by an angle ofK/2) and the integral lines of Q. are parallel 
curves to each other 

The proof of this proposition simply results from introducing the angle so that 
Q.{x,t) = (cos(0(f,?)), sin(0(f,r))) and noting that satisfies the "transport equa- 
tion" 



This property implies that the knowledge of Q. on the cluster boundaries suffices 
to know Q. everywhere inside the clusters. 

The integral curves of Q. provide a mathematical description of the animal files 
in the herd. These curves being parallel to each other, they are consistant with the 
intuition and the observation of animal files in a herd (see fig. [2l). 

The pressure p satisfies an elliptic equation. Indeed, by taking the divergence 
of the equation (12.201) and after easy computations, we get 



where Tr is the trace of a matrix and the exponent T denotes the transpose opera- 
tor. This equation can be equivalently written: 



- {Q^ ■ V,^)2p - (V.r ■ V^)p = -TT{{V^a){V^af), (2.24) 



and only involves the operator (Q ■ V^) applied to p. Since the integral lines of 
Q.-^ are straight lines, equation (12.241 ) is just a one-dimensional elliptic problem 



^;,2 0-(tan0)5.„0 =0. 



Vriiid-a®Q.) v^p) = Tr( ( v.^ti) {y^af) , 



(2.23) 
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Figure 2: Left: schematic figure of a congested zone, where the arrows design the 
vectors Q.. Right: picture of a sheep herd. 



for p posed on this straight line. Knowing the boundary values of p where this 
straight line meets the boundary of the cluster allows to compute p everywhere on 
this lines and consequently inside the cluster (see fig. [2l). Hence, once Q. is known 
inside the cluster, the resolution of this equation only requires the knowledge of 
the boundary conditions for p at the boundaries of the cluster. 

To close the system, i.e. to determine how the solution in the congested do- 
main evolves, we need to determine these boundary conditions. They are not given 
by the formal limit and, in order to determine them, we need to explore another 
route. For this pupose we look at the solutions of the Riemann problem for the 
perturbed and limit systems. Note that if we abandon the constraint of constant 
norm \Q.\ = 1, the non conservative term {Q.0Q.)Vj;p in the momentum conser- 
vation equation (12.131 ) drops out, and we recover a conservative model expressing 
mass and momentum equation. Then, the Rankine-Hugoniot conditions across 
the boundary between the compressible and incompressible regions provide the 
boundary conditions for the pressure at the cluster boundary. The constant norm 
constraint prevents from using this strategy. Therefore, we need to find a different 
route to specify these boundary conditions. 

2.4 Conditions at the boundary of the clusters 

To find the boundary conditions on the cluster boundaries, we need to extract more 
information from the perturbation system (I2.8I) - (I2.10I) than the mere limit system 
(I2.12l) - (|2.15l) . As such, this system is underdetermined. The strategy is to extract 
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Figure 3: Notations at the interface. 



such information by passing to the limit e ^ in some special solutions of this 
system. To underline the difficulty resulting from the non-conservativity, let us 
first look at the Rankine-Hugoniot conditions. We have the following proposition. 

Proposition 2.2 1. If p and Q. are smooth on both sides of a dicontinuity line 
r, then we have 

[p(t2-H-a)]r = 0, (2.25) 
where n is the normal to T and o is the speed ofT. 

l.IfQ. is smooth (i.e. across T and p is smooth on both sides ofT, then 
we have 

[p]r(f2-H^)=0, (2.26) 
where n^ is a unit vector tangent to T. 

The proof of this proposition is omitted. The second relation provides us informa- 
tion when the mean velocity is not tangent to the cluster. In this condition, if the 
mean velocity is continuous, the pressure is also continuous. This implies that the 
pressure is zero on a cluster boundary if the mean velocity is continuous. This fact 
will be supported by the forthcoming analysis. However, as regards the interface 
dynamics, such an analysis is incomplete because the second equation supposes 
that Q. is continuous. 

So as to capture the correct boundary conditions for the pressure p and the 
velocity f2 at a cluster boundary where p and Q. may be discontinuous, we con- 
sider a one dimensional problem in the normal direction n to this boundary (cf. 
figure O. In order to justify this simplication, we introduce the coordinate system 
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(xi.xj) in the normal and tangent direction to the boundary. The angle 6 is de- 
fined so that Q.{x,t) = (cos(0(x,f)),sin(0(x,f))) in this basis. System (|Z81)-(|Z91) 
then becomes (the index £ is omitted): 



(9fp + (p cos 0) +^;,2(P sine) =0, 

[dtO + (cos0(9.^,0 + sin0(9.„0) + (-sin0(9.^j£p(p) +cos0«9.,-2£p(p))] 



- sin0 
cos 6 



) 



(2.27) 



We suppose that all quantities have locally smooth variations in the direction tan- 
gent to the boundary and we focus on the possible sharp variations or discontinu- 
ities in the normal direction. To analyze this situation, we perform a coordinate 
dilation in the xi direction and in time: = 8xi , = 5x2, t' = dt, with 5 ^ 1 . In 
these new variables, all x\ and t derivatives are multiplied by 1/5. Letting 5 — 0, 
we are led to the following one-dimensional system with xi = x: 



Hyperbolic systems like (I2.29I) - (I2.30I) have analytical solutions which are those of 
the Riemann problem. These solutions are associated to initial conditions which 
consist of a discontinuity between two constant states. We will construct the so- 
lutions of the Riemann problem for system (|2.29I) - (I2.30I) and analyze their limits 
as e — > 0. This analysis will give rise to jump conditions at the cluster boundaries 
for these solutions. We will then postulate that these jump conditions are generic 
and valid for all solutions of the limit problem (|2.12| )- (|2.15l ). 

As underlined above, the non-conservative form of system (|2.29I )- (I2.30I) in- 
duces a lack of information about the jump conditions across a boundary. In order 
to waive the ambiguity, we have to make further assumptions. One of them is to 
consider the following conservative system as a way to select discontinuities 



where ^(cos 0) = — In | tan(0/2)| and ^»(cos 0) = — In | sin0|. It is the simplest 
conservation form that system (|2.29I )- (I2.30I) can take. It is obtained by dividing 
(12.301) by sin^ 0. The functions ^ and 4> satisfy: 



dtp + dj,{pcose) =0, 
(9f0+cos0(9;f0 + £sin 



i^ed,,£p{p)=0. 



(2.29) 
(2.30) 



dtp + d_x{p cos 6) =0, 

<9,^(cos(0)) + (9.,($(cos 0) + ep{p)) = 0. 



(2.31) 
(2.32) 




sin^ ' 



1 




COS0 



(2.33) 
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Other conservative forms of (|2.29l )- (|2.30| ) do exist (see appendix |B]) but we con- 
sider this form because it is the simplest. Note that this conservative form is not 
equivalent to the original form (|2.29I )- (I2.30I) because ^(cos0) is an even func- 
tion of 6 . Hence it does not provide information on the sign of . However, this 
information will easily be recovered at the end. We remind that, if all conser- 
vative forms are equivalent for smooth solutions, they differ for weak solutions. 
Therefore, the choice of a particular conservative form must be made on physical 
considerations. Such physical considerations are not available here. In front of 
this lack of information, the choice of the simplest of these conservative forms 
seems to be the most natural one. 

Classical hyperbolic system theory will enable us to solve the Riemann prob- 
lem for (12.3 1| )- (|2. 321) and to take the limit e ^ of these solutions. The limit 
solutions will satisfy some jump relations which we will assume generic of all 
solutions of the limit problem (I2.12I) - (I2.15I) . We now present the result of this 
analysis for such generic solutions. We call "unclustered" region (UC) the do- 
mains where < p < p*, by contrast to vacuum (V) where p = or clusters (C) 
where p = p*. 

Formal Statement 1 The boundary conditions at cluster boundaries or vacuum 
boundaries of system d2.i2D - d2J31) are as follows: 

• Interface (C)-(UC). The pressure jump is given by 

, , [^(cos0)] [pcos0] 

= ^ \^ ~ t^^'^''' ^^-^^^ 

and the shock speed is given by the Rankine-Hugoniot relation 

a = [pcos(0)]/[p], (2.35) 

where the angle brackets denote the jumps across the interface. We note 
that puc = <^fT^d that specifying [p] actually specifies the boundary value 
of p at the cluster boundary. 

• Interface (UC)-(V). The interface speed o is equal to the fluid normal 
speed o = cos = Q.-nat the boundary of the (UC) region 

(J = (cos0)t/c, (2.36) 
and the pressure p is identically zero. 
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• Interface (C)-(V). The interface speed is equal to the normal speed cos d — 
Q.-n at the cluster boundary and the boundary value of p is zero 

a = (cos0)c, pc = 0. (2.37) 

• Interface (UC)-(UC). This is a contact discontinuity between two regions 
of different p. The normal velocity is continuous and equal to the speed of 
the discontinuity 

[cos0]=O, a = cos0, (2.38) 
and the pressure is identically zero. 

We note that all these statements are consistent with proposition 12.21 Section |3] 
provides the detailed analysis which leads to these relations. The dynamics of 
the interface between two clusters (C)-(C) does not follow from the analysis of 
the Riemann problem. We provide a separate analysis of it by introducing the 
so-called cluster dynamics. 

We note that (C)-(UC) interfaces or contact discontinuities (UC)-(UC) may 
incorporate a flip of the sign of in the velocity jump. This has no influence on 
the boundary values of p at the cluster boundary which is the quantity we wish to 
determine by this analysis. The Formal Statement [U is illustrated in figure ID 

2.5 Clusters dynamics 

We now focus on the interface (C)-(C), i.e. a collision of two clusters. The proce- 
dure using limits £ ^ of the Riemann problem does not lead to any conclusion 
since the pressure becomes infinite. Note that this is also the case when dealing 
with the same limit in the standard Euler problem. Therefore, we have to find 
another strategy than using the Riemann problem. We turn our attention to the 
collision between two clusters of finite size and we show that the pressure in- 
volves a Dirac delta at the time of the collision. Such an analysis is inspired by 
the sticky block solutions presented in [10]. 

Consider two one-dimensional clusters which collide at a time t^ (see fig. [5]). 
Before collision, the left (resp. right) cluster at time t < t^ extends between a^{t) 
and b(^{t) (resp. ar{t) and br{t)) and moves with speed 

cos di = a'dt) = b'({t) (resp. cos 0, = a^{t) = b'^{t)), (2.39) 

After the collision, the two clusters agregate and form a new cluster at time t > tc 
extending between a{t) and b{t) and moving with speed cos0 = a'{t) = b'{t). 
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p 



dvc Pvc — 



COS ft, 



(a) Interface (C)-(UC). pc and s are given 
by (l2J4l )- (l235l) . 



P 



(b) Interface (UC)-(V) 



Pc=0 cos Or. 



P 

P' 



COS 0,, 



0n 



(c) Interface (C)-(V) (d) Interface (UC)-(UC) 

Figure 4: Interfaces 
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Therfore, p and are given for t < tc by 

P = P*^at{t)Mt)]^ P*har{t)Mt)V ^ = ^e-Mai{t)Mt)] + ^'"1^(0^(01' 
and for t > tc by 

P = P*l[a(r),MO]' ^ ^ ^l[a(?),MO]- 

where 1/ denotes the indicator function of the interval / (i.e. = 1 if x G / and 
otherwise). We denote by m = b({tc) = ar{tc) the collision point. We look for 
a pressure written as p{x,t) = n{x)5{t — tc). The following proposition provides 
conditions for such type of solutions to exist. 

Proposition 2.3 1- Supposing that p{x,t) = 7t{x)d{t — tc) where 7t is continuous 
and zero outside the clusters, then and K satisfy 

(^(cos e ) - ^(cos e^) ) (m - a{tc) ) + (^(cos 0) - ^(cos 0,) ) {b{tc) - m) = 0, 

(2.40) 

( (^(cos0)-^(cos0f))(7M-Jc) 
7j:{x) = < +(^(cos0)-^(cos0^))(^(/'c)-m), ifxe[a{tc-),m\, (2.41) 
[ (^(cos0)-^(cos0,))(Z7(?c)-^), ifxe[mMtc)l 

2 - Under conditions A2.40\l - iS.41\) . {p,9,p) is a solution (in the distributional 
sense) o/' d2J7]) -f l232D . 

The proof of this proposition is developed in appendix O 

2.6 Conclusion of the analysis 

The underdetermined problem (|2.12| )- (|2.15l ) must be complemented with the For- 
mal Statement [H which determines the boundary values of p at cluster boundaries 
and by proposition [2]3] which determines the evolution of two clusters when they 
meet. Strictly speaking, proposition 12.31 only gives the collision dynamics of two 
clusters in dimension 1 . In dimension 2, clusters may have complicated shapes. 
So, the collision dynamics of two clusters in dimension 2 is a complicated problem 
which will be examined in a future work. At the present stage, problem (12.121) - 
(12.151) complemented with statement [T] fully determines the dynamics of the limit 
system as long as two clusters do not meet. 

A rigorous theory of the well-posedness of system (|2.12l) - (|2.15l) complemented 
with statement[T]is outside the scope of the present paper. Let us just mention how 
a time discretized version of the problem can be computed. Suppose that p"{x). 
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i^"(x),p"{x) are approximations of p{x,t'^),Q{x,t"), p{x,t") at time = nAr. We 
solve the implicit system 

^^^^ + Vf (p"+it2"+i) =0, (2.42) 
_ a« ^ ^ + (Id - ® V3e^"+i = 0,(2.43) 



A? 



with 



r2"+^ = ii±Ai_^. (2.44) 



This form guarantees that p = p = 1 (by taking the dot product by 
and using that \Q."^2 1 = 1). is determined by solving the elliptic equation 

-Vr{{U-a"+-2®a"+^-)V^p"+^)=V:^-{{a"-V^)a"+^^). (2.45) 

on every connected component of the cluster region defined at time t hy {x G 
M? |p''+^(x,f) = p*}. This equation must be supplemented with suitable bound- 
ary conditions on p at the boundary of the cluster. These boundary conditions 
are actually given by the Formal Statement [H with right-hand sides evaluated at 
time t"^^. The resolution of this equation guarantees that V ■ = on every 
connected component of a cluster, and shows that p"+^ = p" = p* on such a clus- 
ter. Of course, the implicitness of the discretization leads to a nonlinear stationary 
problem, and the question of the existence of solutions for such a problem is not 
clear. However, intuitively, it seems that the prescription of the boundary values 
of p at cluster boundaries through the Formal Statement [T] leads to a well-posed 
problem, at least as long as two clusters do not meet. 



3 The one-dimensional Riemann Problem 

3.1 Methodology 

To find out jump relations satisfied by the solutions of the system (|2.12| )- (|2.15l ), 
the strategy is to solve the Riemann problem of the one-dimensional perturbation 
system (12.3 1I) - (I2. 321) and to take the limit e ^ of its solutions. This strategy was 
successfully adopted for a model of traffic jams in [6]. 
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We note that the eigenvalues and eigenvectors of the hyperbolic system (12.311 )- 
(12321) are 



A^(p,0)=cos0±v/ep'(p)p|sine|, F^(p,e)=( ) . (3.46) 



In this conservative system, the domain of is restricted to the interval ]0, ;r[. But 
this is not a problem since our main concern is to find the missing conditions on 
p at the cluster boundary, and these only depend on jump conditions as functions 

of COS0. 



3.2 Solutions to the Riemann problem for (I2.31I) - (I2.32I) . 
3.2.1 Genuinely nonlinear fields 

The Lax theorem provides the local entropic solutions of the Riemann problem 
provided that all the fields are totally genuinely nonlinear (VA^ ■ 7^ 0) or totally 
linearly degenerate (VA|. -r^ =0). Unfortunately, the following result implies 
that the fields are genuinely nonlinear except on a one-dimensional manifold. 

Proposition 3.1 1. The linearly degenerate set (VA|. ■ = Oj consists of two 
curves 'tf^ (each of them corresponds to one characteristic field): 

^l = {{p,e), pe[0,p*l cotane = TG'{p)}. 

where 

r^rnV- ^ (p"(P)P + 3/(p))P ^ iP*-py-^ 



2. For 7=1, the linearly degenerate set tends to the straight lines {0 = 0} 
and {0 = 7t} as £ tends to 0. For 7 > 1, ( resp. j is a one to one, onto 
mapping from [0,p*] to [;r/2,;r] (resp. [0,7t/2]) for all e, called 0/j±(p)- 
For a fixed 6 G]0,2;r[, the inverse map Pid±{Q) satisfies: P* — pf^ ±{d) — 

O(er^). 

The proof of this proposition is easy and is omitted. Thus, the Lax theorem is 
valid at least locally in the neighbourhood of all the states except those which are 
on the one-dimensional manifolds. The second part of the previous proposition 
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shows that all the states have locally genuinely nonlinear fields as e tends to 0. 
Indeed, even if the state converges to a congested state as e ^ 0, its convergence 
is like 0(£^/^) (cf. (12.1 11) ). which is slower than the convergence of the linearly 
degenerate field when 7 > 1 . Therefore, there exists e' such that for all e < e' the 
fields of the converging state are genuinely non-linear. It is also trivially the case 
when 7 equals 1. According to standard nonlinear conservation theory [40], (for e 
small enough) the solutions of the Riemann problem consist of two simple waves 
(shock waves and/or rarefaction waves) of the first and second characteristic fields, 
separated by constant states. 

3.2.2 Shock and rarefaction waves. 

A shock wave between two constant states (p^, 9() and (p^, 9^) travelling with a 
constant speed o satisfies the Rankine-Hugoniot relations: 

[pcos(0)] = a[p], (3.47) 
[4>(cos(0)) + £p(p)] = a[^(cos(0))], (3.48) 

where [/] '■= fr — fe. denotes the difference between the right value and the left 
value of any quantity /. By eliminating o in these equations, we get a non-linear 
relation between the left and right states: 

He{pe,ei,pr,er) := mcos{e))+ep{p)][p]-mcos{e))]\pcos{e)]^0. 

(3.49) 

With a fixed left state, the zero set of He is called the Hugoniot locus and repre- 
sents all the admissible right states, connected to this left state by a shock wave. 

Proposition 3.2 The Hugoniot locus consists of two Hugoniot curves asso- 
ciated to the two car acteristic fields. 

1. The Hugoniot curve M'l associated to Xl_ (resp. to X^) is strictly in- 
creasing ( resp. strictly decreasing ) in the (p , ) -plane. Let h^_ : ] {h^_ ) ^ ^ (0) , :;r [ 
[0,p*[andh^_^ : ]0, (/?*}_) ^ ^ (0) [0,p*[ be the Hugoniot curves asfunctions 
ofG on their domains of definition. 

2. The Hugoniot locus tends to the union of the straight lines {9 = Of } and 
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The proof of this proposition is developed in appendix iDl 

A rarefaction wave is a continuous self-similar solution {p{j):9{j)). It sat- 
isfies the diffential equation 

p'W \ _ F^(p(5),r](5)) 



V'{s) J VA^(p(.),T](.))-r^(p(.),r](5))' 

where rj = ^(cos(0)) is the conservative unknown. Therefore, (p,T]) belong to 
the integral curve of r^. By changing the parametrization of the integral curve, 
we obtain 

p' = ±p|sin0|, 0' = -a/£p'(p)p| sin0|, 
and then the following integral equation 



Q-Q^ = ^ [\ ^^^du. (3.50) 

Jpi V u 

It defines two integral curves issued from the state {p£,0(). The following 
proposition summarizes their main properties. 

Proposition 3.3 1. The integral curve ofr^_ (resp. ofr\) is strictly in- 
creasing ( resp. stricly decreasing ) in the {p,6) —plane. Let i^_ : ] (/l ) ^ ^ (0) , ;r [ 
[0,p*[ and i\ :]0, {i%)^^ (0)[~^ [0,p*[ the rarefaction curves as functions 
ofO on their domains of definition. 

2. For all 7 > 1, the rarefaction curves tend to the union of the straight lines 
{9 = di} and {p = p*}. Moreover, for e]0£,7t[ (resp. d g]0, de[), p* - 

i^_{e) = 0{£r^) (resp. p* - i%{e) = 0{£r^)). 

3. Suppose that the state pf is such that p| p* and ep{pf) — > p(. For all 
p < pf, {i±)^^{p) satisfies: 

\{i%)-\p)-er\<\{i%)-\0)-er\ = O{e^y). 



The proof is developed in appendix |El 

Entropy conditions. In order to satisfy the Lax entropy condition, each 
Hugoniot curve is restricted to right states which have a smaller associated 
eigenvalue than the left state. 
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Proposition 3.4 The eigenvalue (resp. X^) is a decreasing function of p on 
the Hugoniot curve J^f^ ( resp. an increasing function of p on J^_^ ) for < 

cotan^^ (^{-\/y/£p'{p)p^ (resp. 6 > cotan^^ (^\/ y/ep'{p)p^ ). 

Proof Let g : p G [0,p*] g{p) E [0, n] be an arbitrary function. The variation 
of on the graph of g is given by 

VA^ ^=±{x'np)smd+g'{p)i-sme±x'{p)cose). 

where X*^(p) = ^/£p'{p)p. Since x'{p) is positive and the Hugoniot curve {h'L)~^ 
is increasing, is a decreasing function of p on this curve for 6 g]0, :;r[ such 
that (sinO +x^{p)cos0) > 0. Similarly, since the Hugoniot curve (h^)^^ is de- 
creasing, is a increasing function of p on this curve for 6 E]0,n[ such that 
(-sine + ;t'^(P) cos 0) > 0. ■ 



So, in the limit £ ^ 0, the reachable right states are those belonging to the up- 
per half-domain. We denote by = Jt^ n {(p, 0), A|(p, 6) < A^(p^, 0^)} the 
shock curves. 

Concerning the integral curves the admissibility conditions select the 
curves with increasing eigenvalues and so the curves on the lower half-space. 
Therefore, the rarefaction curves satisfy i?^ C ^± n {(p, 0), g]0, 0f[}. The 
union of the shock and the rarefaction curves form the forward wave curve 
Sj^URj^, while the union of their complementary sets form the backward wave 
curve Wt' = ^l\S\ U 

3.2.3 Solutions to the Riemann problem 

Given a left state (p^, 0^) and a right state (p^, 0^), an entropic solution is found 
by intersecting the forward 1-wave curve Wi' issued from the left state and the 

backward 2-wave curve W*'^ issued from the right state (cf. fig. [6l). In the fol- 
lowing study, the curves indexed by - (resp. by -i-) are implicitly those issued from 
the left state (resp. from the right state). Because of the monotony of the shock 
and rarefaction curves, we can classify the different solutions according to the 
positions of the left and right states in the (p, 0) -plane. The following theorem 
describes the solution of the Riemann problem for small e > and is illustrated 
in figure 111 
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Figure 6: Wave curves wl'^ for the left state (p^, Oi) = (0.8, 7t/2) and W^'^ for the 
right state (p^, 0,-) = (0.6,2;r/3). In dashed green lines: the rarefaction curves. In 
continuous red lines: the shock curves. In dotted black lines: linearly degenerate 
sets, p* = 1, 7= 2. 
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Theorem 3.5 Considering a left state (pi, di) and a right state {pr, 0,-), and for £ 
small enough, the solution is given by one of the four following cases: 

1. Case 0( = dr. If p£ < pr [resp. pf > pr], the solution consists of a 1 -shock 
[resp. 1 -rarefaction] connecting {pi.di) to (p,0) (with p G ]p^,pr[ and 
6 > 6{ = 6, [resp. p e ]Pr,P^[ and G < Gi = G,]) and then a 2-rarefaction 
[resp. 2-shock] connecting (p,0) to (p^, 0,). This is summarized in the 
following diagram: 

I „ . shock ,~ ~. rarefaction , _ , . . 

{plM) ^ (P,0) ^ (Pr,0r) lfPl<Pr 

rarefaction ^ ~ shock ^ n\ ^ 
(pe^Ge) — ^ (p,0) — > {pr,Gr) ifpi>pr 

2. Case Gi > Gr (cosG( < cos Gr). The solution consists of a 1 -rarefaction 
connecting [pi.Gi) to (0,0) (with p < pe^pr and G G ]Gr,Gi[) and then a 
2-rarefaction wave connecting (0, G) to (p^, Gy). We get the following dia- 
gram: 

Q\ rarefaction ~ vacuum ~. rarefaction 
{pi,Ge) — > (0,0) — ^ (0,0) — > {pr,Gr). 

3. Case Ge < Gr (cos Gc > cos Gr). There are two sub-cases: 

• if Pr < {h^)^^{Gr) and pf < {h^)^^{Gf), the solution consists of a 
1-shock connecting (p^, Gi) to {p,G) (with p > pi.pr and G ]0^, G([) 
and then a 2-shock connecting (p^G) to (pr, Gr). The diagram is: 

, „ X shock ,~ ~. shock , „ x 

(P^,0^) > {p,G) > {Pr,Gr) 

• if pf > {h^)^^{Gf) [resp. pf > (h^)^^ (Gf)], the solution consists 
of a 1-shock [resp. 1 -rarefaction] connecting {pi,Gi) to (p,0) (with 
p E]p£,pr[ and G > Gr [resp. p G]pr,p^[ and < Gg]) and then a 2- 
raref action [resp. 2-shock] connecting (p, 0) to (pr, Gr). The diagram 
is as follows: 

I „ X shock i~ ~. rarefaction , ^ 
(pe^Gf) — > (p,0) — > {pr,Gr) ifpe<pr 

rarefaction ~ shock , ^ \ -r 

(pl.Gi) > (p,0) > ipr.Gr) lfpl>Pr 
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The detailed proof of this theorem is developed in appendix |F.1[ Let us provide 
some ideas of the proof. For finite e, there exist four kinds of solutions depending 
on what parts of the curves w['^ and W^'^ meet. So, for a fixed left state, the 
state-space is divided in four subdomains. These subdomains depends on the left 
state. However, reminding that the limit of the Hugoniot and integral curves are 
straight lines d = diox p = p* (cf. propositions l3.2| and |3.3l ) for all left states, the 
four subdomains have the same behaviour as £ ^ whatever the left state is. 

3.2.4 The sign of 

The conservative system (12.3 1I) - (I2. 321) does not determine the sign of 9 (if Q is 
supposed to be in ] — k,k\). As mentioned above, this is not important since our 
main goal is to provide connection conditions on p between the left and right 
states. However, it is desirable to determine it in the present analysis, for the 
sake of completeness. For this goal, we cannot use (12.321) because ^(cos0) is 
an even function of . Again, we are facing an indetermination due to the non- 
conservative character of the system. One possible solution is to introduce a con- 
tact discontinuity from 9 io —9 with propagation speed cos0 in the domains 
where p is constant and cos 9 is continuous. If we add such a contact wave, there 
is only one possible construction given by the following: 

Proposition 3.6 Suppose that 9(, 9r G [— TT, k] and 9(, 9r have dijferent signs. 

1. In the subcases cos 9^ = cos 9r and cos 9f > cos 9r of theorem 13.51 the only 
one possible contact wave in the domains of constant p and continuous 
COS0 is located inside the intermediate state and the propogation speed 
equals cos 9. 

2. In the subcase cos 9i < cos 9r, the possible contact waves are those located 
in the vacuum domain. There is no uniqueness of the propagation speed but 
since this contact discontinuity occurs in the vacuum p = region, we may 
consider that 9 is not defined in this region. 

The proof of this proposition can be found in appendix |F.2[ Two cases of the 
Riemann problem with 0^ < < 0^ are represented in Fig. [8l Note that the position 
of the contact wave does not depend on £. So their limits as e goes to zero are 
easily obtained. 
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(e) Case 9c < Or, second subcase, pc < Pr (f) Case 0^ < Or, second subcase, pt > Pr 



Figure 7: Solutions to the Riemann problem for small £ > 0. 
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t 



0e]0^,|0,|[A^ 



cos d 



(a) Case 0^ > 



(b) Case 0^ < |0,.|, first subcase 



Figure 8: Some solutions to the Riemann problem for small £ > and —K< Or < 



3.3 The solutions of the Riemann problem in the limit e ^ 

In order to study the limit £ ^ 0, we introduce converging sequences of left and 
right states 



and we look for the limits of the solutions of the associated Riemann problems. 
There are three cases to consider: either none of the two states tends to the con- 
gested state (pe.pr < p*), or one of the two does (p^ < p*.,pf p*) or both of 
them do (p|,p^ p*). The case (p| P*:pr < P*) is obtained by symmetry 
from the case (p( < p*.,pf. p*): the left and right quantities have to be ex- 
changed and the arrows have to be flipped (like in the first case of theorem [331) . 
Since the solutions of the Riemann problem are bounded and monotonous, all the 
sequences belong to a bounded subset of 5V(]R) and consequently, to a compact 
subset of L|jj^(M). So we only need to prove the uniqueness of the limit of con- 
verging sequences to prove the convergence of the whole sequence and we can 
consider that the convergence is in the almost everywhere sense (up to the extrac- 
tion of a subsequence). 

As a guideline, we mention that, compared with the system with finite £, the 
limit Riemann problem has two additional properties: the appearance of clusters 
which corresponds to the saturation of the constraint p < p* and the disappear- 
ance of rarefaction waves and their transformations into contact waves. In the 
subsequent statements, the term "limit" is a short-hand for "limit of the solution 
to the Riemann problem of (12.3 1D - (I2. 321) " as £ ^ 0. 



((P|,0|),(P^0.')) 




((P£,0,),(P,,0,)) 
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3.3. 1 Case p£ < p*,pr < p* (see fig. ID) 

Proposition 3.7 (Case p£ < p* , pr < p*) There are only three cases: 

(a) Subcase 0( = G,. The limit consists of only one contact wave connecting 
{pi,ee) to {pnOe): 

, „ X contact f „ X 

The travelling speed is equal to cos Bi. 

(b) Subcase Of > Or. The limit consists of two contact waves connecting the two 
states to a vacuum state: 

, n N contact „ X Vacuum „ x contact , „ x 

{pl.Oi) ^ (0,0^) > (0,0,) ^ (pnOr). 

The travelling speeds are respectively equal to cos Oi and cos Or- 

(c) Subcase Oi < Oy. The limit consists of two shocks connecting the left state 
{Pi, Oi) to a congested state (p*, 0,p) and then connecting (p*, 0,p) to the 
right state (p,-, Or): 

(P,,0,) (p*,0,p) iPr.Or). 

where 6 is the unique solution of 

E [min(0^,0,),max(0^,0,)], 
and p is given by 

_ ^ mco.WMpco.(e)], i^^^^^j^j ^ ['P(co.(a))] [pco.(a)]. _ ^^^^^^^^^^^^ 

[P\e [P\r 
The shock speeds are given by the Rankine-Hugoniot condition for the density 

Note that in clustered region, since p — p*, the state is determined by the values 
of and p. This is why we add a third component giving the value of p to the 
vector defining the state in the clustered region. 
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(a) Cas 0£ = e,- (b) Cas 9^ > Br (c) Cas 0^ < dr 

Figure 9: Limit solutions of the Riemann problem for e = and pg^pr < p*. 



In this proposition, the quantities [f]£ '■= f — f(, [f]r'-=f — fr denote the dif- 
ference between the intermediate value and the left (or right) value of the quantity 
/ and [f]\'-= fr — fe denotes the difference between the right and left values of 
the quantity /. 

This proposition covers several kinds of interfaces described in Formal State- 
ment [U the case (a) is an occurence of an interface (UC)-(UC), the case (b) of 
an interface (UC)-(V) and the case (c) of an interface (C)-(UC). Moreover, the 
proposition implies that 6 is continuous inside (UC) domains. 

Note that all the intermediate states are explicitly given or are solutions of a 
non-linear equation and so, are explicitly computable. The proof of this proposi- 
tion is given in appendix lG.il 

3.3.2 Case pe< p\pr = p* (see fig. [H) 

This case is typical of the situation at a cluster boundary. In this case, the main 
new feature is the appearance of declustering waves as limits of the rarefaction 
waves. These declustering waves are instantaneous cancellations of the pressure. 
The following lemma details this statement: 

Lemma 3.8 (Limit of rarefaction waves, declustering wave) Let (p^, Or) be a se- 
quence of right states such that pf ^ p* and £p{pf) ^ pr > 0. Introduce a 
converging sequence of states {p^,G^) lying on the rarefaction curves issued from 
the right states, such that p^ < pf. 

If p = limp^ < p*, then the rarefaction wave tends to the combination of a 
contact wave between the state {p,0r) and {p*j6r,Pr) with speed cosOr = A+ 
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and a declustering wave, i.e. a contact wave with infinite speed which cancels the 
pressure, which provides a transition between (p*, Qr^Pr) and (p*, 0r,O). 

Ifp= lim = p* then the rarefaction wave tends to a shock wave with infinite 
speed between the states (p*, Or^p) and (p*, Or^pr), where p = limejc»(p^). 

The proof of this lemma is developed in appendix lG.2[ 

The next proposition provides the solutions of the limit Riemann problem and 
figure [To! schematically describes them. 

Proposition 3.9 (Case p£ < p*,pr = p*) There are only three cases: 

(a) Subcase di = dr. The limit solution consists of one contact wave connecting 
the left state {p£, Gi) to an intermediate congested state (p*, 0r,p = 0) and 
then a cluster contact (with infinite speed): 

(P,, 0,) (p*, 0,,O) ^^^-l""'- ip\ dr^P). 

(b) Subcase 6( > Or. The limit solution consists of one contact wave connect- 
ing the left state to vacuum, and then another contact wave connecting the 
vacuum to a congested and pressureless state (p*, dr,0) and finally a cluster 
contact connecting (p*,0^,O) to {p*,Or,p): 

I „ X contact „ X vacuum „ x contact / * ^ declust. r * n -\ 
{p£,Oe) — ^ (0,0^) — > (0,0,) — > (p ,0„O) (p ,0r,p). 

(c) Subcase Gi < Or. The limit solution consists of one Shockwave connecting the 
left state (p^, O^) to an intermediate congested state (p*, cmd one con- 
tact wave with infinite propagation speed connecting (p*, Or,p) to the right 
state (p*, 6r,p): 

(Pi, de) {p\ ee,p) (p*, er,p), 



where the intermediate pressure p is equal to 

p = [^(cos 0)], _ [4>(cos 0)], . 



The shock speed is given by the Rankine-Hugoniot condition for the density 
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e,,p*,o 



cos 9f 




(a) Case 0^ = Or 



(b) Case B( > Or 



i)nP%pXl 

(c) Case 0( < Or 



* n£ 



Figure 10: Limit solutions of the Riemann problem for e = and p£ < p*, p 

p\ 



In practice, when instantaneous waves occur (i.e. with infinite propagation speed), 
it means that the initial data of the Riemann problem does not spontaneously ap- 
pear during the dynamical evolution of the limit problem. They have to be ignored. 

Like proposition 13. 7[ this new proposition covers several kinds of interfaces 
described in the Formal Statement [T] cases (a) and (c) are occurences of inter- 
faces (C)-(UC), the left wave of case (b) is an occurence of an interface (UC)-(V) 
whereas the right wave of the case (b) is an interface (C)-(V). 

The proof of this proposition is in appendix IG. 31 

3.3.3 Case pi = pr = p*-,pf < pf (see fig. El) 

We assume in addition that ep{pf) and ep{pf.) have finite positive limits, denoted 
hy p(> and pr > 0. Figure [H] provides a sketch of the solutions. 

Proposition 3.10 (Case pi = pr = p* .pf < pf) There are only three cases: 

(a) Subcase Gi = Gr. The limit solution consists of a uniform constant state 
{p\Gi,p(). 

(b) Subcase Gi > Gy. The limit solution consists of two contact waves and two 
cluster contact with infinite travelling speed: 

/ ^* n - \ declust. / i, „ „s contact, „ „^ vacuum, „ contact , „ declust. , „ _ n 
(P ^Ql^Pl) ^ (P ,0^,0) > (O,0£) > (0,0;-) > (p ,0,.,O) > (p ,Gr,Pr), 
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(c) Subcase Of < Or. The limit solution consists of two shock waves with infinite 
propagation speed connecting the left state {p*,Oi,p(') to (p*,0,+oo) and 
then (p*,0,+oo) to (p*, 

{p\0,,pe) {p\ 0, +00) ^ (p*, 0r,pr), 

where is the only solution of 

[^(cos(0))],[cos(e)]. _ /pA-r 

[^(COS(0))],[COS(0)], \prj ■ ^ ■ ^ 



These solutions display only one kind of interface among those discussed in the 
Fomial Statement[il the case (b) is an occurence of an interface (C)-(V). Accord- 
ing to the cases (a) and (c), the solution inside clusters is continuous. However 
the case (c) does not provide a meaningful solution since the pressure becomes 
infinite and this is why Formal Statement[T]does not allow to decide what happens 
at the interface (C)-(C), i.e. a collision of two clusters. It seems to result from the 
fact that in this case the Riemann problem models the collision of two infinite one- 
dimensional clusters. Section [231 provides a description of the collision between 
finite-size one-dimensional clusters. We have seen that the pressure p involves a 
Dirac delta in time. Indeed, according to case (c), the infinite propagation speed 
of the waves inside clusters implies the discontinuity of the function in time: 
= Of + {0 — 0()H{t — tc), where H denotes here the Heaviside function. Then, 
equation (12.321) leads to 

{'¥{d)-^'{Oe))d{t-t,) = d,p, (3.52) 

which justifies to look for a pressure with a dirac delta in time. The Riemann 
problem does not allow to take into account such a pressure. 
The proof of the proposition is deferred appendix IG.4I 

3.4 Connecting the Riemann problem analysis to the Formal 
Statement [H 

We remind that, by contrast with the finite £ system (I2.8I) - (I2.10I) . which is a stan- 
dard hyperbolic system, the limit system (I2.12I) - (I2.15I) exhibits two additional 
characteristics: 
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COS Of " cos Or 

Vacuum 




er,p\pi ' er,p\ph Bt,p\pf^ Qr.p\ph Oe,p*,P( 
(a) Case 9i = Or (b) Case 0£ > 0, 



er.p\ph 

(c) Case 0£ < 9, 



Figure 11: Limit solutions of the Riemann problem for £ = and pf.pf p*. 



Pf < Pr- 



(i) the appearance of clusters corresponding to the saturation of the constraint 
P = P% 

(ii) the appearance of vacuum. 

The previous study helps in understanding the dynamics of the interfaces be- 
tween unclustered states (UC) < p < p*, vacuum (V) p = and clusters (C) 
p = p*. Up to now, no rigorous theory for the limit £ ^ exists and so, we 
cannot have access to these dynamics rigorously. Our method is to investigate 
these dynamics through the inspection of the limit £ ^ of the solutions of the 
Riemann problem of the finite £ system. 

The first remark is that waves with infinite speed correspond to an instanta- 
neous transition from the initial data to some different solution. To some extent, 
this means that the corresponding initial datum is unstable, and therefore, that 
it will never appear spontaneously in the course of the evolution of the system. 
Therefore, we can discard initial data which exhibit this phenomenon, and replace 
them by the one which is found after the infinite speed wave has been applied. 

The second remark is that the various solutions of the Riemann problem can 
be grouped by situations corresponding to the four cases listed in the formal state- 
ment[Il i.e. interfaces (C)-(UC), (UC)-(V), (C)-(V) and (UC)-(UC). As discussed 
in the lines following statement 13. 101 the case (C)-(C) is not accessible by the Rie- 
mann problem analysis because the collision dynamics of two clusters depend on 
their size, and because the Riemann problem only allows to consider infinite size 
clusters. This is why cluster collisions are analyzed separately in proposition 12.31 
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However, cluster collisions is a complex phenomenon in 2D and proposition 12.31 
only provides a one-dimensional analysis. The two-dimensional analysis is still 
in progress. 

To highlight the link between the Riemann problem analysis and the Formal 
Statement [U we point out which solutions of the Riemann problem correspond to 
which case in the Formal Statement [1} 

Interface (C)-(UC): it appears in 



- prop. 13.71 subcase (c): we see that the intermediate congested state is sep- 
arated from the left and right states by two (C)-(UC) interfaces. We notice that 
(12.341) and (12.351) are respectively the relations for the pressure and the shock speed 
stated at prop. 13.71 

- prop. 13. 9[ subcase (c) (if ignoring the contact wave with infinite speed). We 
can also view subcase (a) as a particular case where the velocities of (C) and (UC) 
are equal (in which case = in the cluster and the interface moves with the 
common velocity). Again, these two cases are consistent with (|2.34|) and (|2.35l) . 

Interface (UC)-(V): it appears in 

- prop. 13. 7[ subcase (b) as two contact waves between the unclustered left and 
right states and the vacuum intermediate state, 

- prop. 13. 9[ subcase (b) where the first wave is a contact between the left 
unclustered state and the vacuum middle state. 

In both cases, the velocity of the interface is that of the non- vacuum states and, 
of course, the pressure is identically zero. Therefore, the situation is as depicted 
in Formal Statement [T] 

Interface (C)-(V): it appears in 

- prop. 13. 9[ subcase (b), where the second wave is a contact wave between the 
vacuum middle state and the right clustered state. We notice that in this case, the 
clustered state must have zero pressure (otherwise, a declustering wave instanta- 
neously relaxes the pressure to zero), 

- prop. I3.10[ subcase (c), where the left and right clustered states are seperated 
by a vacuum intermediate state. Again, in this case, the pressure inside the clusters 
is identically zero. 

In both cases, the velocity of the (C)-(V) interface is that of the cluster. There- 
fore, the situation is as depicted in the Formal Statement [T] 
Interface (UC)-(UC): it appears in 

- prop. I3.7[ subcase (a). We see that this situation is that of a standard con- 
tact discontinuity for the uncongested system. The velocities on the uncongested 
states are equal and equal to that of the interface, and of course, the pressure is 
identically zero. Therefore, the situation is again as depicted is the Formal State- 
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ment[TJ 



We feel that these observation provide a very strong support to the Formal 
Statement [T] As pointed out above, this statement allows to close the system 
(I2.12I) - (I2.15I) at least until clusters meet. In the one-dimensional framework, 
propo sition 12 ■ 3 1 provides the cluster collision dynamics. The investigation of clus- 
ter dynamics in the two-dimensional case is still work in progress. 

4 Conclusion 

In this paper, we have studied a continuum model describing a particle system 
with short-range repulsive and long-range attractive interaction. This is a model 
for the study of gregariousness among mammal species for instance. We have 
focused on the effect of the short-range repulsion and looked at the regime where 
the interaction is turned on suddenly when the local density becomes close to 
some limit associated to congestion. We have modeled this effect by introduc- 
ing a perturbation parameter £ and studied the limit e ^ 0. We have shown that, 
in the limit regime, the congested regions are domains where the flow is incom- 
pressible. The complete determination of the limit system requires the knowledge 
of the interface conditions at the boundaries of the congested regions. We have 
derived these conditions by looking at a model one-dimensional situation (corre- 
sponding to the normal direction to the interface) and analyzing the solutions of 
the Riemann problem for the perturbation model (with finite £). Taking the limit 
e — > in the solutions of the Riemann problem allowed us to provide the missing 
conditions at the interfaces. 

The perspectives of this work are, at the theoretical level, to try to provide more 
solid justifications to these interface conditions and to analyze the cluster collision 
dynamics, which is not accessible by the Riemann problem. At the numerical 
level, we will seek numerical methods to solve this constrained hyperbolic prob- 
lem and we will perform numerical comparisons between the continuum model 
and the more fundamental particle system. 
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A Derivation of a macroscopic model of short-range 
repulsive and long-range attractive interactions 

A.l Individual Based Model with speed and congestion con- 
straints. 

We consider A'^ particles in M? labeled hy k E {1, .., A'^}. These particles are discs 
of radii d. The motion of the particles is described by the time evolution of their 
positions and velocity vectors (0^. Like in the Vicsek algorithm [16,21,44], 
the velocity magnitude of each particle is the same, is constant in time and sup- 
posed equal to c > 0. The velocity direction cbk belongs to the unity circle = 

|«eMM«p = l|. This is a usual assumption in the modeling of several bi- 
ological systems like flocks of birds [3], schools of fish [22, 24] or herds of 
sheep [36,37]. 

We start with a simple continuous-in-time model of a particle system subject to 
attractive-repulsive binary interactions which describe the aggregation of particles 
with occupation constraints. The evolution of the positions and velocities is given 
by: 

= ccbk, (A.53) 



dt 
d(Ok 
dt 



(Id v,«^,«-v[^,0, (A.54) 



where v^^" and v[^'' are the attractive and repulsive forces respectively. The 
matrix (Id — cbk® cOk) is the orthogonal projector onto the plane orthogonal to cbk. 
It is applied to both forces in order to keep the magnitude of the speed constant 
in time. and are the local centers of mass of the particle distribution inside 
interaction discs centered at with radii respectively equal to Ra and Rr 

la _ 3\^rM<Ra ir 3\Xrh\<R,- 



I 1 ' El' 

j,\Xj-Xu\<Ra j]Xj-Xu\<R, 

and v[ are scaling factors which provide the intensities of the forces. The 
repulsive force radius Rr is supposed much smaller than the attractive force radius 
Ra- The resulting force attracts the particles towards the center of mass of the 
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particle distribution at large distances and repels them from the center of mass of 
the particle distribution at short distances. To some extent, it is an implementation 
of the attractive-repulsive scheme proposed by Couzin [16]. 

Finally, we suppose that is a constant and depends on the local density 
inside the repulsive interaction disc: 

v« _ V v'' - V (n''\ n'' - ^j,\^j-^k\<Rr^ 

Vk — ^ai Vk — \Pkl 1 Pk — 

where Vr is an increasing function. The function v,- prevents the local density 
from exceeding the maximal density p* which corresponds to the case where all 
particles are in contact with their neighbours. Clearly, p* is the ratio of the maxi- 
mal occupied surface in a disk of radius R, by disks of radii d and is of the order 
of unity. Therefore, the function Vr tends to infinity as p[ ^ p*. We defer the 
explicit choice of the function to the end of the section. 



A.2 Mean-field model, hydrodynamic limit and macroscopic 
model 

The goal of this appendix is to provide a model for large systems of interacting 
particles according to (|A.53I) - (|A.54|) at large time and space scales. For this pur- 
pose, we will perform a sequence of rescalings. The first rescaling aims at taking 
into account the large number of interacting particles: it leads to the so-called 
mean field model. Assuming that the system is included in a fixed box, the limit 
N — > -t-oo implies that the area Kd^ occupied by each particle tends to like 1 /N 
in such a way that the total area Nnd^ occupied by the particles remains constant. 
We denote by a = liniN^+ooNnd^ the fraction of the surface occupied by the par- 
ticles. The second rescaling is a hydrodynamic scaling where large time and space 
scales are considered. Both scaling are classical and are well detailed in [13, 19]. 
They have been applied to swarming model in [1 1, 21, 26]. 



A.2.1 Formal derivation of the mean-field model 

We refer to [42] for classical references on the mean-field limit. We consider the 
empirical distribution f^{x,G),t) defined by 

1 ^ 

fix,(b,t) = - £ 5{x-Xi,{t))5{(d,(bk{t)). 

k=i 
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5{x) denote the Dirac delta on M?, while 5{cb, (do) denotes the Dirac delta on 
centered at G)o (i.e. COq) is the probability measure supported by {wo})- It is 
an easy matter to check that satisfies the following kinetic equation 

+ ca>- W,f + . ( - f!') f) = 0, 

where and F^ are the attractive and repulsive forces, given by 

F^{x,G),t) = Vfl(Id- g3® F^{x,a),t) = (Id- w ® G))|^^, 



with 



■^.7 ,^ _ IUy-x){y-x)p^{yj)dy J K,.{y-x){y-x)p'^{yj)dy 



JUy-x)p^{y,t)dy ' ^ ' ' ^ /i^,(3;-f)p^(3;,Orfy 
A^;rJ^/ij:,(3;-x)p^(3;,OJ y^ 
jK,{y-x){y,t)dy 



where p^{x,t) = J^^gi f^{x, ^,t)dQ. is the local density and Ka (resp. A"^) is the 
indicator function of the disc of radius Ra (resp. Rr). Here, it is clear that more 
general kernels Ka, Kr can be used. 

The formal mean-field limit A'^ — +0° of this model is (we recall that a = 

limTv^+ooA^TTJ^): 

drf + ca)-V^f + V^-{{Fa-F,)f)=0, (A.55) 

Fa{x,G),t) = Va{ld- (0(^G))^a, F,.(x, «, = Vr(Id - W ® (A.56) 

j Ka{y-x)p{y,t)dy j K^iy -x)p{y,t)dy 

Rigorous justifications of this limit are outside the scope of this article. 



A.2.2 Hydrodynamic scaling 

In order to select the relevant scales, we first rewrite our system in dimensionless 
variables. We consider a space scale xq (typically the range R, of the repulsive 
force) and we choose a time scale to = xo/c. The associated dimensionless time 
and space variables are t' = t /to and x' =x/xo- We also introduce scaled collision 
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kernels K'^ ^ such that Ka^ri^ox') = K'^^^ix), scaled intensities = Vac/x^, = 
Vrc/x^ and a scaled distribution function /' = af. In the case where Ka^r are 
indicator functions of balls of radii Ra,r, this amounts to rescaling the radii to new 
values R'^ = Raj/xQ. After removing the primes, the system in the new variables 
and with the new unknowns is similar to (IA.55I) - (IA.58I) but with c = 1 and a = 1. 

To derive the large time and space dynamics, we introduce the following 
change of variables x = rjx, i = rjt with rj <^ I. In the new variables, the dis- 
tribution function /"^{x, cbj) = f{x, (b,t) satisfies the following system (omitting 
the tildes): 

r^{dJ^+CD- V,/^ ) + ■ ( (FJ' - F,^ ) ) = 0, 

llKjl('-f){y-x)p^{y,t)dy 



-> rj -> Tj 

Fj'{x,cb,t) = vj]{ld-a)'^a))^r , {x,a),t 



^ !K^{^-f)pn{y,t)dy 

llKj^(fif)p^mdy\ r 

, \ / .X , P'{x,t) = / f'^{x,a>,t)d(b, 



where kJ] and kJ^ are the scaled interaction kernels and , , the scaled inten- 
sities. 

We first suppose that the repulsive kernel KJ^ and the repulsive intensity v,^ 
are unchanged in the scaling: Kr = K,, Vr (p) = v^(p). This means that the range 

of the repulsive force is supposed of order rj . To analyze the limit T] ^ 0, we first 

-> rj 

need an expansion of c,r in terms of rj . The following lemma provides the result 
for an isotropic kernel Kr {Kr{z) = Kr{\z\)) 

Lemma A.l Under suitable regularity assumptions on p^, we have the expansion 



y-x 



p'l{y)dy = ap'^{x)+o{r]), 



„,'i/'f'(¥)p''(?'"''>\ 
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where a = J Kr{\z\)dz andB = J Kr{\z\)z®zdz= {J Kr{r)r^drld). 

The proof of this lemma is elementary and omitted. In the case where Kr is the 
indicator function of the disc of radius R, , the coefficients a and B are equal to 

a = nR^ and B = ;r^Id. Now we consider the scaling of the attractive kernel Ka 
and attractive intensity Vq . We suppose that the attractive force remains non-local 
as 77 tends to and weaker than the repulsive force. To express these assumptions, 
we suppose that the scaled attractive kernel Ka and intensity are given by 

K^{z)=KMz), yJ! = ri\,. 

For simplicity, we choose = 1 . We also fix the space unit xq in such a way that 
B/a = 1. In particular, in the case where is the indicator of the ball of radius 
Rr, we can fix xq = Rr/2. 

Under all these modelling assumptions and thanks to lemma IATT] the system 
can be written formally, in the limit rj ^0: 

S,f+a- V-J + Vs,-((F,:-Fr)f)^<3, (A.59) 

F„(x,s,,) = (Id - ffl0S)f„, ur^,,) = ( x^ifziiiyi^piv 

V }Ka{\y-x\)p{y,t)dy J 
F,(^,«,0 = v,(p)(ld-w®w)|,, t{^,t) = ^j^^. (A.61) 

A.2.3 Macroscopic model 

The last step is to obtain the dynamics of macroscopic quantities associated to 
the flow. Here we will only consider the density and momentum. We find that 
under suitable regularity and decay assumptions on /, the density p = J fd(b and 
momentum pQ. = J fcbdco satisfy the following system of mass and momentum 
balance equations: 

dtp+VrP^ = Q. (A.62) 

V .1 Ka[\y-x\)p{y,t)dy J p{x,t) 
To close system (IA.62I) - (IA.63I) . we assume that / is a monokinetic distribution: 
f{x, cb,t) = p{x,t)8{G),Q.{x,t)), (A.65) 
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with \ Q{x, t)\ = \. This assumptions presupposes that a local equilibrium is reached 
where all particles are locally aligned. Although no justification of this assumption 
can be made at this point, the features displayed by the system seem meaningful 
in view of gregariousness modelling. We find: 

dtp+V^pQ. = 0, (A.66) 

dt{pa)+v^-{pa®a) = p{id-a®a){^a-yr{p)t), (a.67) 

where and are given by (IA.64I) . Factoring out p in (IA.67h . using (IA.66I) . we 
also get the following form of the system: 

dtp+V^-pa = 0, (A.68) 
+ V^Q. + v,.(p) (Id - f2 ® = (Id - f2 ® a)^a, (A.69) 

^ ^ ' V jKai\y-x\)p{y,t)dy J' ' ' p{xj)^ 



A.3 Repulsive force intensity and macroscopic model 

Let us return now to the choice of the function v^. This function tends to infinity 
when p ^ p*. Like in the traffic model devised in [6] , we assume that this function 
behaves like when p <^ p* and tends to infinity when p ^ p*. The prototype 
of such a function is 

P(P)= . ^ .y , (A.71) 

[v~p) 

where 7 > 1 . We will keep this example constantly in the paper for simplicity 
but the results are valid for all functions having the same properties. We consider 
that Vr{p) = pp' (p). In this way, we suppose that repulsion acts like a standard 
presure force in a gas, but, when the density reaches the congestion density p*, 
the pressure tends to infinity. Since the equation for Q. is used instead of that 
for pQ., the interpretation of p in standard gas dynamics terms would rather be 
that of an enthalpy (i.e. p'{p) = P\p)/p where P is the actual fluid mechanical 
pressure), but the results would be similar if we considered the equation for pQ. 
instead. Indeed, because of the constraint \Q.\ = 1, the system is non-conservative 
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in the projection term (Id — f2 ® Q). Finally, we get the following system 



dtp + V^-pQ = 0, (A.72) 

dta+a-v^a+{id-a^a)v^p{p) = (id-a^a)^,,, (a.73) 

V lKa[\y-x\)p{y,t)dy J 

This system provides the starting point of the present article. Since this paper 
is focused on the treatment of congestion phenomena, we remove the non-local 
attractive force. Indeed, this term is a zero-th order derivative term and does not 
intervene in the jump relations across discontinuities. 



B Conservative laws for the one-dimensional system 

In this appendix, we are looking for conservative forms of the one-dimensional 
system (|2.29I )- (I2.30I) . The most general conservative form is written: 

drg{p,e) + dj{p,e)=o. (B.75) 

where g and / are smooth functions of p and 6. The following proposition ex- 
hibits an infinite set of such conservative forms. 

Proposition B.l If{gjf) is a conservative form of ^2.29\iA2.30\i . then their partial 
derivatives are related by 

^ = -^cose--^sin0y(p), .^ = -^cos0--^psin0. (B.76) 

dp dp de de de dp^ 

Moreover, if g is a function with separated variables g{p, 6) = u{9)v{p), then u 
and V satisfy 

pv"{p)=kp\p)v{p), (B.77) 
u"{d) -f {cotand)u{d) = ku{9), (B.78) 

where k is a constant real number Each e M gives rise to possible {g, f) pairs. 

Proof Performing the chain rule in (IB. 751) and using (12. 291) . (12. 301) . we easily get 
(IB.76I ). Then, using that differentiations with respect to p and commute, (IB.76|) 
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gives rise to an elliptic equation satisfied by g and inserting the hypothesis of 
separated variables, we obtain (|B.77| ). (|B.78I ). 

Equation (IB.78I) is the Legendre differential equation (in polar coordinates). 
The two-dimensional vector space of solutions of this equation is spanned by the 
Legendre functions of first and second species and each of them gives rise to 
possible (g,/) pairs. □ 

The solutions of (IB.77I) exists for all G M. However, they have a priori no ex- 
plicit expression exept for k = 0. In this case, the 2-dimensional vector space of 
solutions of (IB.77I) is spanned by {l,p} and for (IB.78L is {l,^(cos0)}. We can 
actually check that the following (g,/) pairs 

U,/) = (p,pcos0), (B.79) 
(gj) = (^(cos0),4>(cos0) +pip)), (B.80) 
(gj) = (pvi'(cos0),pcos0^(cos0) +P{p)), (B.81) 

where P is an antiderivative of pp'{p), are non trivial solutions. The conservative 
form studied in this article corresponds to the pairs ( IB. 791) and (IB. 801) . The pairs 
( IB. 791) and (IB. 811) form another such conservative system. 



C Proof of proposition |23| (cluster collisions) 



Proof 1- Let xq e [a{tc),m] and h G '^^{D"), where D" is a neighbourhood of xq 
in D (cf. figure [5]). We apply the Green formula on the domain D": 

< ^f^(cos0)+(9.^4>(cos0),/7> = - 11 ^{co^e) dxh + ^f {cos e)dthdtdx 

h[{^{cosd),^'{cose))-n]ds 

dD" 

xo+ri 



(^(COS0)-^(COS0^)) / h{tcM'^) 



where < .,. > denotes the duality brackets. Since we look for p{x,t) = n{x)5{t — 
tc) , we also have 



— <dxP,h>= // pd^hdtdx= I 7t{x)d^h{tc,x)dx= — dx7t{x)h{tc,x)dx. 
JJd" JxQ—r\ Jx(^~r] 

(C.83) 
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If (12.321 ) is satisfied, then we have 

< (9f^(cos0) + a,$(cos0),/z >= - < dxp,h >, (C.84) 
and equations C.82I ) and (IC.83D imply that 

(^(cos0) -^(cos(0£))) = -dM^o)- 
The same arguments (for any x G [a{tc),b{tc)]) lead to 



(^(cos0)-^(cos(0^))), ifxG [a{tc),m\ 
(^(cos0)-^(cos(0,))), ifxG [m,b{tc)] 



and (supposing ;r continuous) to 

( (^(cos0)-^(cos(0£)))(m-x) 
;r(x) = < +(^(cos0) -^(cos(0r)))(Z7(?c) -m), ifxG [a(/'c),m], 
[ (^(cos0)-^(cos(0,)))(M?c)-x), ifxG [m,M?c)], 

Supposing that p and then ;r equal zero outsidse the clusters, we get 

;r(a(?,.)) = ('i'(cos0)-^(cos(0^)))(m-a(?,)) + (^(cos0)-^(cos(0,)))(Z7(?,.)-/? 

2 - Let hhe a test function in the neighbourhood D' of D (cf. figure [5]). We 
denote Di = Dr\{t <tc} and D2= Dr\{t >tc}. Applying Green's formula, we 
obtain: 



< dtP + dx{p cos G),h >= — jj p cos Gdxh + pdth dtdx 
h[{p cos 9, p) ■n]ds+ / h[{p cos 9 ,p) ■ n]ds 

JdDo 



dDi JdD2 

fbiitc) 

{-p*cos{9i)+ai{t)p*)h{t,ai{t))dt+ / p*h{tc,x)dx 

tc-8 Jai{tc) 



E 

i<E{i,r} 



+ r {-p*cosi9,)+b\{t)p*)h{t,b,{t))dt 

.hc-8 

ftc+S rb{t,) 

+ / {-p* COS 9 +a {t)p*)h{t,a{t))dt+ / p*h{tc,x)dx 

Jtc Ja(tc) 

- r^\-p* COS 9 +b'{t)p*)h{t,b{t))dt 

Jtc 

—p* / h{tc,x)dx — p* / h{tc,x)dx + p* / h{tc,x)dx 

J aeitr) Jar(tr) J aitr) 



0. 
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since nds = ±{—l,x'{t))dt on the left and right sides of the domains Di and D2 
and nds = ±{0,l)dx on theirs top and bottom sides. The last equality stems from 
the identity 

(beitc) - ai{tc)) + {br{tc) -ar{tc)) = {b{t,) - a{t,)). 

The density equation (12.311) is satisfied in the distributional sense. 

If we now apply Green's formula with a test function h E ^^{D), we obtain 

< dt'i'icose) + dj,^{cose), h >=- JJ $(cos0)a^/z + ^(cos0)(9f/z dtdx 
h[{^{cose),'¥{cose))-n]ds+ [ h[{^{cos 0),'i' (cos G)) ■ n]ds 

dDi JdD2 



fin foyjc) 

/ (^(cos0)-^(cos0|))/i(r,,Jc)Jx+ / (^(cos0) -^(cos0,))/i(r,,x)Jjc 

Ja{tc) Jin 

rb(tc) 

— j dxTth{tc.,x)dx = — < p,h > . 

J ait A 



'a(tc) 

Eq. (|2.32l) is satisfied in the distributional sense. Note that in this case, the test 
function has a compact support in D since ^(0) is not defined in the vacuum 
region p = 0. ■ 



D Proof of proposition |33l (study of the Hugoniot 
loci) 

In this section, we provide a detailed study of the Hugoniot curves. Let (p^ , di) G 
]0,p*[ X ]0, k[ be an arbitrary left state. We need to find the geometric behaviour 
of the Hugoniot loci associated to this left state. The classical theory of nonlinear 
conservation laws provides only information on the local behaviour of M'^. Each 
J^, consists of a one-dimensional manifold tangent to the integral curves 
of the right eigenvectors up to the second order. In the (p, ^(cos 0))-plane, the 1- 
Hugoniot curve is thus locally decreasing and the 2-Hugoniot curve M",^ is lo- 
cally increasing because of the direction of the vectors = (±p | sin 1 , £p'{p)p). 
In the (p,0) -plane, the 1 -Hugoniot curve J^_^ defines a locally increasing func- 
tion 6 = {h^)^^(p) while the 2-Hugoniot curve defines a locally decreas- 
ing function = {h^)^^{p). Actually, this property is global (i.e. (ht)^^ (resp. 
{h^)~^) is a globally increasing (resp. decreasing) function of p for all p g]0, p* [). 
To prove this, let us begin with a simple and useful lemma. 
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Lemma D.l For all u ^ [—1, 1], the function : v g] — 1, 1 $(v) — u^{v) is 
convex and has a minimum at the point u. In particular, we have 

Vv ^ M, (4>(v) - 4>(m)) - u (^(v) - ^(i/)) = Uv) - Uu) > 0. 

The proof is elementary and omitted. We now analyze the behaviour of in 
more detail when e becomes small. Proposition [3]2] is an immediate consequence 
of the following lemma. 

Lemma D.l The behaviour of Jif^, does not depend on the left state. Let 
(Pf , di) be a left state. Then: 

(i) Suppose Gr is fixed. The function py H£{p£,G£,prj6r) has at most two 
zeros and there exists e' > such that for all £ < e', the function Pr 
HE{pg, Gg, pr, Gr) has only one positive zero. This zero tends to p* as £ tends 
toO. 

(ii) Suppose pr is fixed. Then Ve > 0, the function Gr — > H£{pi,G(,pr,Gr) has 
two zeros, one lower and one larger than G(, and both of them tend to G( as 
£ tends to 0. 

The Hugoniot locus tends to the union of the straight lines {G = Gf} and {p = p * }. 

Note that these results imply that the Hugoniot locus consists of two monotonous 
curves as functions of p (otherwise H with fixed Gr would have more than two 
zeros). The local behaviour of the Hugoniot locus enables us to determine that the 
increasing curve is associated to the first eigenvalue and the decreasing curve 
to the second eigenvalue A^. 

Proof (i) Let us fix the left state (p^, Gg) and the right angle Gr. So as to get a 
more readible proof, the function pr H^^pi, G£,pr, Gr) will be denoted by H but 
its derivative will be denoted by a partial derivative dp^ H. We look for the zero set 
of H in the interval ]0, p* [. We compute: 

-J^ipr) = e {2p\Pr) +P"{pr)pr) > 0. 

As p and its first two derivatives are strictly positive on ]0,p*[, the function H 
is strictly convex and thus has at most two zeros. Moreover, the value of H at 
pr = p( is strictly negative, 

//(P^) = -P£[^(cos0)][cos0] <0, 
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if dr is not equal to Og. Like the function p, H tends to +00 when p,- tends to the 
maximal density p* . Then H has only one zero in ]p^, p*[. We have 

H{0) = {ep{pt) + [^(cos 0)] cos 0^ - [4)(cos 0)]) . 

Lemma ID. II implies that the second term of this expression is strictly negative 
and thus H{0) becomes stricly negative for small e. Thanks to its convexity, we 
deduce that there exists e' such that for all e < e' , the function H has no zero in 
the interval ]0,p^[. 

To show that the only zero of H tends to p*, let us rewrite H as follows 

H{p,) = [£p{p)] [p] + [4>(cos 0)] [p] - [^(cos 0)] [p] cos 0, - [^(cos 0)] [cos 0] P£, 

and thanks to lemma jPTl the zero of H satisfies 

[ep(p)] [p] = - [4)(cos 0)] [p] + [^(cos 0)] [p] cos(0,) + [^(cos 0)] [cos 0] pi 
> [^(cos0)][cos0]p^>O. 

So we can easily conclude that the zero of H tends to p*. 

(ii) Like in the first point, let us denote the function Or H£{p£, 0e,Pr, •) by 
H. First, the value taken by H at Or = Of is positive: 

H(0,)=£[p(p)][p]>O. 

Some easy computations leads to the following expression of the first (partial) 
derivative of H: 

dH 1 

3^(0.) = {Pe [cos 0] + [^(cos0)] p,sin2 0,) . 
oor smy,- 

As p£ and pr are positive, the sign of the derivative is the same as the sign of 
[cos 0]. Thus H is increasing on [0, 0^] and decreasing on [0^, n]. Moreover using 
the fact that ^(m) =^{u) + log( 1 + m) , we can write H as 

H{er) = 4>(C0S 0,.) [p (1 - COS 0)] + log(l + COS 0,.) [p COS 0] + A(£, p,., 0,) 

where A is a bounded function. It implies that H tends to —00 when Or tends to 
0. In the same way and by using the identity ^(m) = — ^'(m) +log(l — m), we can 
show that H tends also to —0° when Or tends to n. We deduce that H has exactly 
two zeros. 

Let us remark that 

H,=H,-{\-e)[p{p)][p]. 

This implies that H^\0) = H^\{\ - £)[p{p)][p]) and then that the zeros of He 
tend to Of as e tends to 0. ■ 
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E Proof of proposition |33l (study of the integral curves 



of the right eigenvectors) 



Proof 1 . We easily check that d'{p) = ^^/£p'{p)/p, leading to the result. 
2. For a fixed p, the quantity 



(4)-i(p)=0,TVif r 



converges to Oi as £ goes to 0. For a fixed 0, the quantity ^/e ^ \J~^^-^du^ ~ 

G — 6(1 is a constant. So as £ tends to 0, the integral term has to tend to 
which implies the convergence of = i^{0) to p*. Besides, the function inside 

the integral behaves like O (^y/£{p* — m)^^ j when p p*. This leads to a 
diverging integral for 7 > 1 . Then the integral behaves like O (^^/£ (p * — pj ) ^ ^ j 

and thus we get p*- 4(6) =o(^£^j. 

3. Let £'>Oandp < pf . From the rarefaction curve equation (|3.50|) . {i±)^^{p) 
satisfies 

|(4r P-0.I</ \ -^^du = V^ ^ ^^du, 

Jo y u Jo (p*-m)^ 

Assuming that the limit of £p{pr) is finite, we get p* — p^ = 0(£^). Thus, the 

7+1 

function inside the integral behaves like \^{p* — u)^^ when p ^ p*. This leads 

^/£{p* — pf.) 2~ 

andthus like 0(£2r). ■ 



F Proofs of theorem |331 and proposition |321 (solu- 
tions of the Riemann problem for £ > 0) 



F.l Proof of theorem 3.5 



Proof Let (p^, Q() and (p^, Qr) be left and right states respectively and let us sup- 

f e 

pose that the intersection of the 1-forward wave curve W_ issued from the left 
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b E 

State and the 2-backward wave curve W^' issued from the right state reduces to 
one point (p,0) (in all the proof, the 1-wave will be implicitly relative to the 
left state while the 2-wave curve will be implicitly relative to the right state). 
The solution of the Riemann problem depends on which parts of the wave curves 
meet: for instance, if (p, 6) is the intersection of the 1-shock curve with the 2- 
rarefaction curve, then the solution will be the combination of a shock wave (with 
a speed given by the Rankine-Hugoniot condition (13.471) ) and a rarefaction wave 
separated by an intermediate constant state. To find where the intersection on the 
shock curves is located, the main arguments will be the monotony of the wave 
curves given by propositions 13.21 and 13.31 (independently of the location of the 
states in the (p, 0) -plane) and their convergence speed to their asymptotic limit. 
The wave curves will be considered as functions of 6 in their domain of definition: 



where w_ is an increasing function and w_|_ is a decreasing function. The functions 
h± and i± are respectively defined in propositions l3.2l and l3.3[ Let us examine the 
different cases suggested by the theorem successively: greater or lower or equal 
to 6,-. For the reader's convenience, the corresponding geometric configurations 
of the wave curves are illustrated in figure [T2l 

Case Oc > 0^ (Fig. [HI (a)). From proposition^ (?1)"H0) (resp. {i%)-\0)) 
tends to Of (resp. 0, ) as e goes to zero (and the third point of the same propo- 
sition asserts that it is still the case when pr tends to p*). So, assuming that 
there exists a such that (?'")"^(0) < (/^)~^(0), there exists /3 < a such that 

since the domains of definition of w_ and 

w+ are respectively [(/^)^^(0),:;r[ and ]0, (z^)^^(O)], the only intersection point 

of and is the intersection of the 1 and 2-rarefaction curves at (?'^)^^(0). 
As e decreases, the intersection point disappears since the domains of definition 
are separated. However, the integral curves meet the {p = 0} axis at the states 
(0, /f.(0)), (0, i^{0)) and these states are connected by vacuum. 

Case 6( < 9,-. We suppose that pi is lower than pr. For all e, the increasing 
1-shock curve issued from the left state divides the domain [p^,p*] x [0^,2;r] (to 
which the right state belongs) in two parts: the left domain where the right state is 
on the left side of the 1-Hugoniot curve issued from the left state {pf < h^_{9f)) 
and the right domain where the right state is on the right side of the 1-Hugoniot 
curve issued from the left state {pf > (0^^)). 

- Assume that for all e the right state is on the right side of the 1-shock curve 




for0G[(/l)-i(O),0£], 
ht for e [0f,:/r[, 




for G [0,0,], 
for0e[0„(/^)-i(O)[, 
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e 




(C) Of < e,-, Pr < /!- (d) Oe = Or 



Figure 12: Schematics of the intersections of the wave curves (proof of theorem 
I3.5|) . Only the parts which meet are represented. 
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(Fig. [T2I (b)). We consider — on the interval [(/l)^^(O), (/*j_)^^(0)], which 
is the intersection of the domains of definition of and w^. The function — 
is increasing. We have (wl — w*^)((«l)^^(0)) = — (w+)((z-)^^(0)) < and 
(wi — w^) (Or) = ( Or) — pr > 0. So the only zero of — is in the interval 
[(il)^HO), 0^]. So the intersection point of the two wave curves is the intersection 
of the 1 -shock curve and the 2-rarefaction curve. This corresponds to the second 
subcase of the third case of the theorem. Note that the limit of the 1 -shock curve 
(proposition l3.2l) implies that the limit of this case should be considered only if py 
tends to p*. 

- Assume that for all e the right state is on the left side of the 1 -shock curve 
(Fig. [121 (c)). Since p£ is lower than p,., the left state is also on the left side of 
the 2-Hugoniot curve issued from the right state: pf < pf = h^{6^) < h^{6^). 
We again consider the increasing function — on its domain of definition 
[(/'l)^^(O), (z*^)^^(0)]. This function is negative at Of and positive at Or. So the 
intersection point of the two wave curves is the intersection of the two shock 
curves and p > p^, py. This corresponds to the first subcase of the third case of the 
theorem. 

If p£ is greater than pr, the decreasing 2-shock curve issued from the right state 
divides the domain [Pr,p*] x [0, 0,] and the same arguments as before lead to the 
result. 

Case 0{ = Or (Fig. [I2l (d)). Assume that pr < p£. We again consider the in- 
creasing function w'L—w^. It is positive atO = 0{ and negative for = (z*i ) ^ ^ (0) < 

Oi (since /z+ ((?i)"HO)) > p^ > 0). So it equals zero for a value G < Or. So (^p, 6^ 
is the intersection of the 1 -rarefaction curve and the 2-shock curve, which leads to 
the solution given in the first case of the theorem. The case pr > p( is similar. ■ 

F.2 Proof of proposition 

Proof 1. Consider the domain where {p = p^}. 

If the left state is connected to the intermediate state via a rarefaction wave, 
then this rarefaction fan is contained between the speeds = cos Oi — ^/ep'{p()pe\ sin Og 
and = cos — y^ep'{p)p | sin 1 . Since < cos Oi, the domain {p = p^ } can- 
not contain the contact wave with speed cos 6(. 

If the left state is connected to the intermediate state via a shock wave, then 
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we have cos < cos 6g and pg < p, which yields 

— cosOg-l — (cos 0^ — cos 0) 

Pe-P 
- pfi ~ 
= cos0H ■ — (cos 0— cos 0^) <COS0, 

p-p( 

So the domain {p = pf} cannot contain the contact wave with speed cos Og. 

So in both cases, the domain {p = pi} cannot contain the contact wave with 
speed cos 0^. The same arguments show that the domain {p = p^} cannot contain 
the contact wave with speed cos 0^. 

Finally, we easily check that a contact wave with propagation speed cos can 
occur within the intermediate domain {p = p}. Indeed, if the intermediate state is 
connected to the left state (resp. to the right state) via a rarefaction wave, then we 
have < cos < A:^ and if the intermediate state is connected to the left state 
(resp. to the right state) via a shock wave then < cos (resp. 5+ > cos 0). 

2. Like in the previous point, the contact wave can be located only in the 
intermediate state {p = 0}. But the propagation speed is not unique: it can be all 
the intermediate speeds between the two fans of rarefaction. ■ 



G Proofs of lemma and propositions 13771 1X91 and 



3.101 (limits of solutions of the Riemann problem) 



We recall that the quantities indexed by - (resp. by +) are implicitly those related 
to the left state (resp. the right state). The characteristic speeds related to the 
intermediate state will be denoted by A^. 



G.l Proof of proposition 1377 



Proof (a) Let us suppose that p£ < Pr (the opposite case is similar). According to 
theorem [33] and proposition l3.3l the intermediate state angle 0^ tends to Of (since 
p^ belongs to the interval ]Pf,Pr[ for each £). In addition, it is easy to check that 
the two speeds and tend to cos 0^ as £ tends to 0. So the rarefaction wave 
turns into a contact wave with speed cos Or. As regards the shock wave, its speed 
is given by 



)^COS 0^ — p^COS 0^ _pCOS0^ — COS0£ 

^ = cos Oe + p ^ 
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and tends to cos dg as e tends 0. Finally, the two waves coincide and make a single 
contact wave. 

(b) As in the first case, it is easy to check that the two rarefaction waves (cf. 
theorem [331) turn into contact waves with speeds respectively equal to cos Or and 
cos 

(c) According to proposition 13.21 h^_{p^) tends to 9^ < Or- So, theorem [331 
implies that we are looking for the limit of the intersection point of the two shock 
curves issued from the left and the right states. These intermediate states (p^, 6^) 
are the solutions of the non-linear systems 

He{Pi,ee,p',e')=0, //e(p^0^p„0,)=O, (G.85) 
p > max(p^,p,), e e ]min(0^,0,),max(0^,0^)[. (G.86) 

From proposition 13. 21 for all G ]min(0^, 0^), max(0^, Or) [, the largest zero of the 
function p — > He{P£, O^.p, 0) tends to p* as £ — 0. Indeed, if p^ does not tend to 
p*, 0^ simultaneously tend to Of and to Or (which is different from Of), which is 
absurd. Therefore, p^ increases and tends to p*. Besides, we have the equality 

ep{p^) [p]e = [^(cos 0)]e[p cos 0]e - [^{cos 0)]e[p]e + £p{p^) [p]e, 

which implies that ep{p^) is bounded as £ tends to (because 0^ is bounded too) 
and we deduce that £p{p^) converges to a non-zero value p. Finally, we can easily 
check that the system given in is equivalent to (IG.86i ). ■ 



G.2 Proof of lemma [X8 

Proof Suppose that A+ = lim A:^ is finite. Since £p(p^) pr>0, then £p'{pf) — > 
+00 and consequently A+ The limit rarefaction wave has a fan for speeds 

s belonging to ] A+, +oo[. The 2-rarefaction wave satisfies, for all 5 e]A+, +°°[, 

s = A+(p(5), 0{s)) = cos(0(^)) + ^£p'{p{s))p{s)\ sin(0(5))|. 

1 

So, for a fixed 5, p(5) ^ p* as £ ^ and we have {p* - p{s)) = C>(£r+i ). Thus, 

£p{p{s)) (= 0{£^^^)) ^ as £ ^ 0. So the rarefaction wave tends to the 
combination of a shock wave between the states (p, 0^,0) and (p*,0^,O) with 
speed A+ and a declustering wave. 

If p < p*, then equals and the previous arguments apply. Let us look at 
the case p = p*. If is finite, then in the previous conclusion the shock wave 
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disappears since the two states on both sides of the shock wave are equal. And it 
confirms that p = limep(p^) equals zero. If is infinite like limA:^, then the 
rarefaction wave turns into a shock wave with an infinite speed between the states 
{p\0r,pe) and {p\0r,pr). ■ 



G.3 Proof of proposition 13.91 

Proof (a) We want to apply lemma lB^ Here the intermediate state is the intersec- 
tion of the 2-rarefaction curve and the 1 -shock curve. The intersection state exists 
for all £ by the monotony of the two curves (cf. propositions 13.21 and 13.31) . Let us 
note that there is no reason to have a finite limit of since the intermediate state 
can tend to p*. By a compatcness argument, we can restrict ourselves to prove the 
uniqueness of the limit of convergent solutions. So let us consider several cases: 

Case (i) — > p^. In this case, the 1 -shock disapears and the solution is given 
by lemma [321 

Case (ii) p^ ^ p G ]p^,p*[. In this case, it is easy to check that the 1-shock 
becomes a 1 -contact and the limit of the 2-rarefaction is given by lemma |3^ 
Case (iii) ^ p*. Now let us look at the limit of £p{p^). We have 

Heipe, Oe, p^ e') = [4>(cos 0) + £p{p)]f [p]g - [^(cos 0)]^ [p cos 0]^ = 0. 

Since cos 6 tends to cos Or (cf. third point of proposition l3.3l) . the terms [p cos 0] 
and [p] are bounded and taking the limit £ — 0, we get 

So either [p]^ tends to or ep(p^) tends to 0. Thus £p{p^) tends to 0. Finally 
lemma \3lS\ applies and we conclude that the 2-rarefaction tends to a declustering 
wave. Now let us look at the limit of the shock speed. It is written 

£ P*^COS0*^ — P^COSQ^ _gCOS0*^ — COS0^ 



s — ^- = cos 0| + p 



Pi P'-'Pi 

and so tends to cos 6,-. So the 1-shock tends to a contact discontinuity. 

(b) The limit of the 2-rarefaction is given by lemma 13.81 (with p^ = and 
A+ = cos 0^).The 1 -rarefaction wave turns into a contact wave as before. 

(c) We first consider the case where the intermediate state is the intersection 
point of two shock curves (for all £, h^{0^) > pf): by the monotony of these 
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curves, is larger than and so tends to p* too. Besides, the intermediate state 
(p^, 0^) satisfies 

HeiPe, Of, p^ e') = [4>(cos 0) + £p{p)]^ [p]g - [^(cos 0)]^ [p cos 0]^ = 0, 
p^ >max(p(,,p,), 0^G]0£,0,[, 

which implies that £p{p^) converges to a value denoted p. By taking the limit 
e — > in the Rankine-Hugoniot relation (13.481) . we obtain 

[^(COS 6)1 [p cos e]r = [4)(C0S 0) + £p{p)l [p]„ 

we obtain 

[^(cos0)]^[cos0], = O 

(0 is bounded and [p]r tends to 0), which implies that equals 0^. Finally, we 
have 

[£p{p)]r = mCOS 0)]r^-^^^ - [*(C0S 0)],. 

If the limit [ep(p)]^ is non zero, then the propagation speed is infinite and if it is 
zero, there is no discontinuity. 

Consider now the case where the intermediate state is the intersection of the 
1 -shock curve issued from the left state and the 2-rarefaction curve issued from 
the right state (for all £, h^_{Of) < pf.). From proposition ^. 31 (point 3), the inter- 
mediate angle 0^ tends to 0^. Besides, thanks to proposition l3.2[ the intermediate 
density p^ tends to p*. From (13.491) . £p{p^) converges to a value denoted by p 
which is given by the limit Rankine Hugoniot relation. So the limit intermediate 
state is (p*, 0,-,p). Finally lemma lXSl applies: the rarefaction turns into a shock. ■ 



G.4 Proof of proposition 13.101 

Proof (a) Since the intermediate density p^ is comprised between the left and 
right ones: pf <p^ < pf (cf. theorem [331) and since pf,pf^p*,we also have 
p^ p*. For the intermediate angle 0^, the previous proof shows that it tends 
to Gr- Let us note that the 2-rarefaction wave tends to a contact wave (because 
^f+T^+ ~^ +°°)- Let s be the limit of the 1-shock speed and let us note that s is 
lower than cos 0^. Like in prop. 13. 9[ subcase (c), the intermediate pressure is equal 
to 

= _ [pcos0],m^ 
[Pli 
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Now let us look at the limit of the shock speed. For finite e, the shock speed is 
given by: 



£ p*^cos0^ — p^cos0^ ^gCosO*^ — COS0^ 



s = = cos Oi + p 



1 

1 



Since ep{pf) ^ pr>0, we have e^' = 0{p* -pf) and then = 0{p*-p^). On 

the other hand, from lemma [33] we have {i^)^^{p) — Or = 0{£^) and therefore 
we get 

cos0''-cos0/ = -2sm — ^ j sm — ^ j = C>(e2r). 

Thus, we easily get that tends to cos Q( and then that the pressure p equals p(. 

(b) Here the proof is similar to the case where only one state converges to the 
congested state (see proof of prop. 13.91 ). 

(c) Consider the case where the solution is the limit of two shock waves. By the 
monotony of the shock curves, the intermediate density is larger than the right and 
left ones (cf. theorem [33]) and so it tends to the congested density too. Suppose 
that the intermediate angle is not equal to 0^. As regards the 1 -shock speed, we 
have _ 

£ [PCOS0]£ ~e , ^ (COS0^ — COS0^'' 



^ = r 1 = COS 0'^+pr 

W P'-P£ 

and so the limit 1 -shock speed is -co. Besides, we have 

[ep{p)]e = [^(cos 0)],!^^^ - [$cos 0],, 

[Pk 

which implies that £p{p^) tends to +co. Then we have 

[p COS0] 



[P]> 



[^(cos0)],= [ep(p)],+ [4)(cos0) 



Since the right hand side tends to +°°, the 2-shock speed has to tend to +oo too. 
If the intermediate angle tends to 0/, then it does not tend to 0, and the same 
arguments apply. The quantities 

[epipMph = [^(cos GMp cos 0]^ - [4>(cos 0)]4p]^, 

[ep(p)]r[p]r = [^(COS 0)],[PCOS 0], - [$(C0S 0)],[p], 
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are bounded. The limit of their quotient is 



[£p{p)Up]r [^(COS0)],[COS0], 

[£pipMp]i [^(cos0)]4cos0]/ 
Besides, it is easily checked that 

[£pip)]r[p]r P-Pr P* - Pr 

[£p{p)]e[p]e e-o p-pf e^O p* -p/ 

where the last equivalence results from the fact that (p* — p) =o(e^) since £p{p) - 
+00 and (p* - p£^r) = O(eT), £t = 0{p* -p(^r)- Finally, we have 

P*-Pr ^ f £p{pe) \ ^(PlY 
p*-pe \£p{pr)J ^ \PrJ 

Consider now the limit of a solution consisting of one shock wave and one 
rarefaction wave. From lemma \3M the intermediate angle 6^ tends to Or and 
from the Rankine-Hugoniot relation, we have 

[£P{P)]1 [Ph = [^(cos 0)], [p cos G], - [4>(cos 0)], [p], . 

So, since 6 is bounded and the densities p^ < p^ tend to p*, the right hand side 
tends to a non zero value: p* [^(cos 0)]^ [cos 6]^. Because [p]^ tends to 0, £p{p^) 
has to tend to +0°, which is absurd since £p(p'^) < £p{pf). ■ 
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